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INTRODUCTION 

At high temperatures, the properties of a gas are influenced 
by vibrational excitation, dissociation, electronic excitation, and 
ionization? all sinks of energy which absorb some of the kinetic 
energy of molecular motion and change the thermodynamic and 
transport properties of the gas rather strongly. Moreover, the 
different modes of energy may not relax to equilibrium in the time 
scale of concern, for example the dwell time of the flow about a 
hypersonic vehicle. Usually, however, the collisions will 
establish a Boltzmann like distribution of states within each 
energy mode very quickly; only a few collisions are required to 
randomize the kinetic energy states in a most probable 
distribution, and to a first approximation the same is true of the 
molecular vibrational states if one neglects the distortion of the 
distribution that occurs due to anharmonicity in the upper 
vibrational levels. Although the exchange of kinetic energy is 
minimal in collisions between free electrons and heavy particles, 
the electron-electron collisions quickly establish a Boltzmann 
distribution of electron kinetic energy states. This state of 
pseudo steady nonequilibrium is described by a gas with three 
different temperatures; T, T a , and T v . These temperatures describe 
the Boltzmann distributions of heavy particle kinetic energy, of 
electron kinetic energy, and of vibrational states respectively. 

Experimentally, the transfer of energy between diatomic 
vibrational states and free electron kinetic energy is known to be 
very rapid. For example, this fact is used to establish 
vibrational population inversions in the CO Electric Discharge 
Supersonic Laser 1 . Moreover, the calculation of diatomic 
vibrational excitation by electron collision has been calculated 
using first principles 2 ' 3 , and the relaxation of free electron- 
molecular vibration states towards equilibrium with one another is 
orders of magnitude faster than the relaxation of vibrational 
states by heavy particle collision 4 . Thus a reasonably good 
approximation for high temperature gases, at least in the time 
scale of concern in hypersonic flow problems, obtains when T. and 
T v are equalized®. 

The relaxation of excited electronic states of atoms and 
molecules toward equilibrium is not as well understood. The 
assumption is often used that excited electronic energy, like 
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vibrational energy, is in equilibrium with the free electron 
kinetic energy. Park 6 published some estimates of rate 
coefficients for electron bombardment excitation of N and 0 atoms, 
based on the work of Lotz 7 , which indicate equilibration between 
free electron energy and excited electronic energy of the atoms is 
indeed rapid, comparable with rate coefficients for excitation of 
vibrational states by electron collision. Thus the approximation 
that excited electronic states are described by T. is no doubt a 
reasonably good one. 

However, once equilibration with electron collisions has been 
achieved, these collisions no longer change the distribution; 
additional energy to the vibrational and excited electronic modes 
feeds in from collision with heavy particles. A portion of this 
energy is in turn transmitted to free electron kinetic energy 4 , and 
in fact the electron-vibration and electron-electronic state energy 
transfers are no doubt the principal mechanisms by which electron 
kinetic energy is eventually brought into full equilibrium with the 
heavy particle kinetic energy of the gas. The direct energy 
exchanges that occur in electron-heavy particle collisions are 
negligible in comparison. 

In the case of vibrational excitation there exist some good 
experiments backed by reliable theories which provide quantitative 
vibrational relaxation rates 8 - 15 , which can be used to predict the 
relaxation of hypersonic flow to full equilibrium at moderately 
high temperature. However, the rates for electronic excitation are 
needed to complete the picture at very high temperature. As 
suggested above, the equilibration between free electron kinetic 
energy and the excited electronic states of heavy particles is 
probably rapid. But the rate of electronic excitation by heavy 
particle collisions is required for prediction of electronic energy 
in the temperature range around 5000°K where normally few free 
electrons exist. When the free electrons do appear, the heavy 
particle excitation of electronic states is a channel that competes 
with vibrational excitation as a mechanism for equilibration of 
electron kinetic energy, and at higher temperatures where the 
molecules disappear the heavy particle excitation of electronic 
energy is probably the principal mechanism for this equilibration. 
Accordingly, the purpose of this paper is to develop a semi- 
classical theory of electron state excitation of atoms and 
molecules caused by collisions with heavy particles, and to compare 
the relaxation rates obtained in this process with the relaxation 
achieved by free electron collisions. 


ELECTRONIC STATE RELAXATION BY ELECTRON COLLISION 

Let the free electron kinetic temperature be T. and the 
excited electronic state distribution be described by the 
temperature T*. Assume that collisions between gas particles in 
excited states transfer electronic energy so readily that the 
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distribution becomes Boltzmann like quickly and that the relaxation 
of this process can be treated as instantaneous on the time scale 
of concern. Also assume that free electron collisions are the only 
effective means of introducing additional excited electronic state 
energy. Then the rate of change in electronic energy E* is 


dE* 
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where n e is the number density of free electrons, n A is the number 
density of atoms in state i with electronic energy e w n is the 
total number density of atoms, and k ±j is the rate coefficient for 
excitation from state i to j by free electron collision. Eq.(l) 
accounts for all upward and downward transitions in electronic 
state over all levels up to some maximum level j_ ax . 


The number densities nj are given by the temperature T* 
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where g A is the degeneracy of the i 1 * 1 state and Q(T*) is the 
electronic state partition function of the atom 

0(T*) - Y, 9i e^ i/kT ' ( 3 ) 
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By detailed balancing at equilibrium with the 
temperature T a , the ratio of the rate coefficients is 
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Thus Eq . ( 1 ) can be expressed 
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dt 
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The deviation of electronic excitation energy per atom from 
its equilibrium value at the temperature T. is 
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The product of electron density n e and the relaxation time r for 
the process is n e times Eq.(6) divided by Eq.(5) 
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Park 6 approximates the calculated values for the rate 
coefficients k 13 with the empirical form 
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The values of the coefficients C 1;) and given in Table 2.1 of 
Park's book 6 were used to calculate the n.r products given by 
Eq.(7). The results for N atom relaxation are shown in Fig.(l) and 
for O atom relaxation in Fig. (2). Note that the relaxation is 
insensitive to the initial electronic temperature T*. The 
relaxation rates are expected to decrease monotonic ly as 
temperature increases. For e-N collisions this expectation is 
satisfied except for very low values of temperature, where the rate 
coefficients given by Eq.(8) are no doubt unreliable. However for 
the case of e-0 collisions, the results look very strange indeed. 
The relaxation rate increases at temperatures above 15000 degK and 
plateaus unrealistically at still higher temperature. Probably 
this is an artifact which results because so many of the excitation 
rate coefficients are essentially zero in Park's table for e-0 
collisions, which seems physically unlikely. Some modification of 
the constants C n not specified in Ref. (6) has been applied by Park 
in Ref. (13), but that is unlikely to compensate for the missing 
rate coefficients in the relaxation rate calculation. However that 
may be, the density relaxation time products are very short, the 
order of lCT 17 to 10* 1 * mol sec/cc. This is even smaller than the 
density relaxation times for electron collision production of 
vibrational states 4 (where n.r is the order of icr 15 to 10 -1S mol 
sec/cc), so rapid equilibration of excited electronic states and 
free electron kinetic energy can be anticipated. The principal 
concern here is how these electron collision rates compare with 
rates for electronic state excitation by heavy particle collisions. 


ELECTRONIC STATE EXCITATION BY HEAVY PARTICLE COLLISION 

The atoms or molecules in high energy collision will be 
assumed to experience a simple exponential repulsion potential 
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where A is a constant which establishes the size of the collision 
potential and L is a scale factor which determines the steepness of 
the potential as the distance R between the collision partners 
changes. Since the major contributions to the transition 
probability will occur where the potential is the steepest, long 
range attractive parts of the potential can be accounted for by 
increasing the collision energy of the pair by the depth of the 
attractive well and adjusting A and L to duplicate the steepest 
part of the potential, which is in effect traversed with the 
velocity of this increased collision energy. The long range 
attractive forces between molecules will be rather negligible at 
high temperatures in any case. However, for the atomic collisions 
a deep attractive well occurs when the electron spins of the outer 
electrons of the colliding atoms happen to pair; then the higher 
effective collision energies will need to be taken into account. 

For the potential of Eq.(9), the solution for the 
intermolecular distance as a function of time is analytic in the 
case of head on collision 11 . 
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where u is the collision velocity in center of mass coordinates, E 
is the collision energy (/iu 2 /2), and R„ is the distance of closest 
approach 
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If one of the collision partners makes a transition from state 
a to state b during the collision, the excited electron produces a 
dipole which creates an electric field V at the position of the 
other particle 


V 


e r cos8 
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where r is the displacement of the electron charge e from its 
nuclear center and 0 is the polar angle the vector r makes with the 


5 


interroolecular distance vector R as shown on Figure 3. A small 
perturbation energy (-aV 2 ) occurs between the particles, due to the 
dipole induced in the inert collision partner (that is, the partner 
not making a transition) , where a is the polarizability of that 
partner. This ripple perturbation added to the interaction of 
Eq . ( 9 ) gives the total interaction for collision when transition 
occurs . 

U - A e ~ R/L - a V 2 (9.1) 


The small perturbation transition probability for promotion of 
state a to b is 
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where <j> a and <p b are the spatially dependent parts of the wave 
functions for the states a and b respectively. The time dependent 
parts of the total wave function provide the circular frequency for 
the transition 




(14) 


The quantum average of the first term in Eg. (9.1) vanishes due 
to the orthogonality of the wave functions, and the resulting 
transition probability can be expressed as a product of several 
dimensionless factors 

p ab ~ ((«/*!) (e 2 /a 0 ft<x>) ( ajR 0 ) 6 H ab F ab f (13.1) 


The dimensionless transition matrix element H ab is 

H ab - < (p; | (r/a 0 ) 2 cos 2 0 | <p a > 

and the dimensionless Fourier transform is 
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POLARIZABILITY OF ATOMS AND MOLECULES 

The polarizability of multielectron atoms and molecules needed 
in Eq.(13.1) is a function of frequency 
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where the sum extends over all optical transitions of the inert 
collision partner, f A is the oscillator strength of transition i, 
©it is the radian frequency of transition i, and y i- s the damping 
coefficient 
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The oscillator strengths needed in Eq.(15) are tabulated for 
atoms by Weise, Smith, and Glennon 14 , for example. 


TRANSITION MATRIX ELEMENTS 

Evaluation of the transition matrix elements of Eq.(13.2) 
requires knowledge of the wave functions <p a and <p b . Although these 
functions are not known precisely for atoms in general, the upper 
excited states at least are nearly hydrogen like, in which case the 
Slater type screened wave functions are a fair approximation. For 
the state of quantum number n 

<P„ - iWzr/n) 3^,(0,*) ( 16 ) 


where R^,, is the usual hydrogenic radial wave function, Y lm is the 
spherical harmonic angular function, and the effective screened 
charge Z is related to the observed energy E n of the n th quantum 
state 



Where the ground state is involved in the transition a 
somewhat better approximation may be required for the ground state 
wave function, such as the partial quantum number functions of 
Bates and Damgaard 15 . Although transition matrix elements 
determined with these approximate wave functions will not be exact, 
the correction factor should be the order of unity so the results 
will at least provide reasonable estimates of the magnitude of the 
transition probabilities and exhibit correct functional relations. 
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FOURIER TRANSFORMS 


The Fourier transform given by Eq.(13.3) and Eq.(10.2) is 
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The subscripts a and b for the specific transition in question are 
omitted for economy in notation. 


Transform the variables to 
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Then the dimensionless Fourier transform becomes 

F-tf e^dz 
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The denominator of the integrand in Eq.(19) has sixfold roots 
at ±z 0/ where 

z D - i arccos(e' a ) (19.1) 


The integral of Eq.(19) is obtained by integrating around a 
half circle in the complex plane, along the real axis from -<» to +<» 
and then back along the half circle at infinity in the upper plane; 
this encloses the pole at +z Q . Then the integral along the real 
axis is just 2?ri times the residue at +z a . 

P f dz 2lI p_£ Res(z 0 ) (20) 

— (1 — In sech z) 6 
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The residue is obtained by multiplying the integrand with 
(z-zj% taking the 5 th derivative with respect to z, and then 
evaluating that derivative at z 0 . This is a somewhat lengthy but 
straightforward process. The result is an analytic function of 
collision energy which can then be integrated over a Maxwell 
Boltzmann distribution of collision energies by numerical 
quadrature . 
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( 21 ) 


Define the functions G(z) and g(z) 
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The fifth derivative is now evaluated at the pole z Q . 
According to Eq.(21) the values of g and its derivatives at the 
pole are 


( 22 ) 
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where the subscript o means that the quantities are evaluated at 
z 0 . The values of g and its derivatives required are thus 
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For simplification, the following positive real factors of 
these derivatives are used 
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In terms of these positive real quantities, the fifth derivative of 
G evaluated at the pole z c is 
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Then the Fourier transform is just the product of 2irRi and the 
residue at the pole 
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EXCITATION RATE COEFFICIENTS 


With the transition probabilities P ab calculated from 
Eq.(13.1), the rate coefficients can be determined- Just as done 
for vibrational transitions, the head on collisions are assumed to 
establish the functional form for the transition probability and 
the absolute magnitude of the rate coefficient is obtained by 
multiplying the result with a constant cross section S of which 
should be the order of 1CT 15 cm 2 for the ground state atoms. For the 
excited states this constant is increased by the range of the 
hydrogen like outer electron wave function 


(l-ejl) 2 


(27) 


where e a is the electronic energy of state a and I is the 
ionization energy of the ground state of the atom. The rate 
coefficient is obtained by numerical integration over a Boltzmann 
distribution of collision energies 

OQ 

k ab - uS 0 f (S a /S 0 ) P ab x e x dx (28) 


where u is the mean molecular velocity in center of mass 
coordinates, x is the dimensionless collision energy in units of 
kT, and x c is the activation energy in these units. S a /S Q and P ab 
are calculated as functions of the dimensionless collision energy 
and the integral is performed numerically by quadrature. 


EXAMPLE RESULT FOR AN N-N COLLISION 

The relations developed above are now used to calculate the 
probability of transition of an N atom in the (3s 2 P) state to the 
(4s 2 P) state when it is perturbed by a collision with an N atom in 
the ( 4 S°) ground state, which is taken to be unchanged in the 
collision. 

The polarizability of the ground state atom is given by 
Eq .(15). The 4 S° state has two strong absorption lines at 83336 cm' 1 


12 


and at 88135 cm" 1 with f numbers of .35 and .13 respectively 14 . 
These lines are so far from the transition frequency <i) ab 

0 ) . - ok-m - 104200-86200 - lBOOOcnT 1 (29.1) 

av d a 


that the damping factor can be neglected. Then 

. e2 Y £i — - 3.386 

al 4 n 2 almc 2 7 


(29.2) 


and the dimensionless factor required in Eq.(13.1) is 
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The matrix element H ab is expected to be the order of 10. This 
factor is approximated using Slater type screened wave functions 
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The effective nuclear charges are adjusted to agree with the 
observed energy levels 

Z 2 - n 2 (I-E) / (e 2 /2a 0 ) 

Z a - [9 (117345-86200) /109571] 1/2 - 1.605 (31) 

Z b - [16 (117345-104200) /109571] 1/2 - 1.385 

To the accuracy of this approximation the transition matrix 
element is 

- Jcos 2 0 sin0 dQ J (R b r 2 R a ) r 2 dr (32.1) 

O O 
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(32.2) 
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Finally the transition probability is a function of collision 
energy, which for this particular transition is given by 
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To evaluate the Fourier transform of Eq.(13.3) the constants 
of the repulsive potential of Eq.(9) must be specified. For the 
purely repulsive potential that occurs most frequently between 
colliding atoms, the constants A = 305 ev and L = .35 A have been 
chosen. These constants fit the potential given by Gilmore 16 fairly 
well; they are assumed to also approximate the repulsion with the 
excited atom except that the size of the cross section is increased 
as in Eq .(27). The transition probability of Eq.(34) calculated 
for this potential interaction is shown on Fig. (4) as a function of 
collision energy. Except at collision energies the order of 30 ev 
or larger these probabilities are small compared with unity so that 
the small perturbation method is justified at low temperatures. 
However, at high temperatures the tail of the Boltzmann 
distribution contributes much to the rate coefficient where the 
small perturbation method breaks down. To compensate for this the 
transition probabilities are limited to unity at collision energies 
above 30 ev. 

The rate coefficient is obtained by integrating the transition 
probability over a Boltzmann distribution of collision energies for 
a range of different temperatures. Fig. (5) shows as a function 
of temperature between 5000 and 50000°K. These rate coefficients 
are considerably smaller than the free electron excitation rates as 
expected. The relation for electron collision excitation of the 
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same transition, that is the excitation of N(3s 2 P) to N(4s 2 P), is 
shown by the dashed line on Fig. (5); Park's expression 6 for this 
rate coefficient is 

k 5 _ 8 - 4.3 • 10 16 | 


The rate coefficients calculated for ground state N atom collision 
excitation of this transition are given as a function of 
temperature in the following table. 

TABLE I. RATE COEFFICIENT FOR N(3s 2 P) + N(*S <> ) - N(4s 2 P) + N(*S°) 


T , degK 

k,cc/mol 

-sec 

log 1Q k 

5000 

1.960 

X 

10' 2 

-1.708 

10000 

1.066 

X 

10 5 

5.028 

15000 

1.708 

X 

10* 

8.232 

20000 

1.472 

X 

10 10 

10.168 

25000 

2 . 688 

X 

10 11 

11.430 

30000 

1.940 

X 

10 12 

12.288 

35000 

7.897 

X 

10 12 

12.898 

40000 

2.354 

X 

10 13 

13.372 

45000 

5.544 

X 

10 13 

13.744 

50000 

1.096 

X 

10“ 

14.040 


This sample calculation illustrates how other transition 
reaction rates could be determined by this method. Although the 
transition matrix elements are not highly accurate, these are not 
the crucial factors that the Fourier transforms are as far as 
establishing the correct magnitude of the transition probabilities. 
Once the rate coefficients have been determined for all the 
important transitions, the relaxation rates can be found from 
Eq.(7) using the Master Equation from Eq.(5) to obtain the rate of 
change in electronic state energy caused by atom-atom collisions. 


CONCLUDING REMARKS 

The relaxation times for electronic excitation due to electron 
bombardment of atoms has been found to be quite short , so that 
electron kinetic temperature T. and the electron excitation 
temperature T‘ should equilibrate quickly whenever electrons are 
present. However, once equilibration has been achieved, further 
energy to the excited electronic states and to the kinetic energy 
of free electrons must be fed in by collisions with heavy particles 
that cause vibrational and electronic state transitions. 


10000 V e -25422/r e cc/mol-sec (35) 

T 
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The rate coefficients for excitation of electronic states 
produced by heavy particle collision have not been well knowh. 
However, a relatively simple semi-classical theory has been 
developed here which is analytic up to the final integration over 
a Boltzmann distribution of collision energies 7 this integral can 
then be evaluated numerically by quadrature. Once the rate 
coefficients have been determined, the relaxation of electronic 
excitation energy can be evaluated and compared with the relaxation 
rates of vibrational excitation. Then the relative importance of 
these two factors, electronic excitation and vibrational excitation 
by heavy particle collision, on the transfer of energy to free 
electron motion, can be assessed. 
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FIG.(l) RELAXATION OF N ATOM ELECTRONIC ENERGY BY 

ELECTRON IMPACT 


FIG. (2) RELAXATION OF O ATOM ELECTRONIC ENERGY BY 

ELECTRON IMPACT 


FIG. (3) EXPONENTIAL COLLISION POTENTIAL PLUS INDUCED 

DIPOLE INTERACTION 


FIG. (4) PROBABILITY OF ATOMIC N(3s 2 P) TRANSITION TO 

ATOMIC N(4s 2 P) DUE TO COLLISION WITH ATOMIC N( 4 S°) 
AS A FUNCTION OF COLLISION ENERGY 


FIG. (5) RATE COEFFICIENTS AS A FUNCTION OF TEMPERATURE 
FOR ATOMIC N(3s 2 p) TRANSITION TO ATOMIC N(4s 2 P) 
DUE TO COLLISIONS WITH ELECTRONS AND ( 4 S°)N ATOMS 
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FIG. (3) EXPONENTIAL COLLISION POTENTIAL PLUS INDUCED 

DIPOLE INTERACTION 






